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Nitrogen-Vacancy (NV) color centers in diamond have emerged as promising quantum solid- 
state systems, with applications ranging from quantum information processing to magnetic sensing. 

One of the most useful properties of NVs is the ability to read their ground-state spin projection 
optically at room temperature. This work provides a theoretical analysis of Purcell enhanced NV 
optical coupling, through which we find optimal parameters for maximal Signal to Noise Ratio 
(SNR) of the optical spin-state readout. We conclude that a combined increase in spontaneous 
emission (through Purcell enhancement) and in optical excitation could significantly increase the 
readout SNR. 


Isolated solid-state quantum systems, such as quan¬ 
tum dots, NV centers in diamond and other point-like 
defects, coupled to photonic degrees-of-freedom, could 
serve as single-photon sources and quantum nodes BB- 
This hybrid approach to quantum networks and quantum 
information processing architectures combines the bene¬ 
fits of solid-state storage and manipulation with photonic 
communications and distribution. 

A commonly encountered problem is the limited cou¬ 
pling efficiency between the quantum system and light, 
which hinders the scalability of the proposed quantum 
networks and the rate of information transfer. Several 
approaches have been suggested to enahnce the coupling 
efficiency including optical cavities and waveguides [SE], 
as well as nanofabricated plasmonic and dielectric struc¬ 
tures BE]. Generally speaking, optical antennas selec¬ 
tively enhance the optical coupling to a specific resonant 
mode, resulting in modified spontaneous emission, called 
Purcell Factor (PF), and in high directionality of this 
emission, leading to a better photon collection efficiency 
m. The resulting enhancement could have a major 
impact on the applicability of these systems, advancing 
toward the goal of efficient solid-state/light interfaces and 
potentially achieving single-shot readout of the quantum 
state of the system. 

However, in certain cases the modified behavior of the 
quantum system in the presence of the resonant antenna 
(and given the resulting Purcell enhancement of its spon¬ 
taneous emission) might degrade its usefulness as a quan¬ 
tum emitter/node. Specifically, for the case of NV cen¬ 
ters in diamond, the optical readout of the defect’s spin 
state relies on the spin-dependent branching ratio be¬ 
tween radiative and non-radiative decay paths. By en¬ 
hancing spontaneous emission, optical antennas may di¬ 
minish the difference in fluorescence emitted by different 
spin states and hence reduce readout contrast. 

In this work we consider an NV center coupled to an 
antenna [Fig. [lja)] and investigate its effect on spin- 
state readout. A leading candidate for the realization of 
such an antenna is a plasmonic structure, since it can 


accommodate a wide resonance to cover a large part of 
the broad NV phonon-sideband emission, yet still allow 
for large Purcell factors (PFs) due to its very small op¬ 
tical mode volume BE]- We focus on the interplay be¬ 
tween the readout Signal-to-Noise Ratio (SNR) and the 
Purcell enhancement of radiative decay for this system, 
and find experimentally feasible optimal parameters that 
significantly enhance both optical coupling and readout 
fidelity. We note that this result relies on the assump¬ 
tion that spin-mixing terms are nonradiative in nature 
mm, since otherwise the improvement in SNR will be 
limited (as further discussed below). 

The NV center consists of a substitutional nitrogen 
atom and a vacancy occupying adjacent lattice sites in 
the diamond crystal. The electronic ground state is a spin 
triplet, in which the m s = 0 and ra s = ±1 sublevels expe¬ 
rience a ^ 2.87 GHz zero-field splitting [Fig.JTJb)], while 
a static magnetic field can further split the ±1 sublevels 
to create an effective two-level system. The NV spin can 
be initialized with optical excitation, detected via state- 
dependent fluorescence intensity, and coherently manip¬ 
ulated using microwaves m 

The optical readout of the NV center relies on exci¬ 
tation and fluorescence intensity measurements. Experi¬ 
mentally, the spin state of the NV is extracted from the 
emitted fluorescence intensity under optical (green) exci¬ 
tation at a wavelength of 532 nm. The excitation laser 
is switched on for spin-state readout, and simultaneously 
fluorescence (in a red phonon-sideband between 650 — 800 
nm) is collected during a readout time window T. The 
dynamics of the NV under optical excitation are deter¬ 
mined by rate equations, based on the transition rates 
between different energy levels as depicted in Fig. [ljb). 

The ability to accurately measure the state of an NV 
center relies on the difference in the number of photons 
emitted during measurements of the two spin projections 
in the ground state, \g,m s = 0) (ground state with spin 
projection m s = 0) and \g,m s = ±1) (ground state 
with spin projection m s = ±1), and is quantified by 
the measurement SNR defined below. Let 7V 0 and N± 
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FIG. 1. NV coupled to an optical antenna, energy levels 
and allowed transition, (a) A general schematic example of 
a plasmonic antenna designed to yield a Purcell enhance¬ 
ment of the radiative rate and to direct the emission, (b) 
Energy levels and allowed transitions for NV center with 
non-radiative spin-mixing transitions in the excited state. 
K e is the spin-preserving excitation rate (with wavelength 
= 532 nm), Kf is the spin-preserving radiative decay rate, 
K s is the non-radiative decay rate to the singlet and Ko is 
the non-radiative decay rate from the singlet. K m is the 
non-radiative non-spin-preserving rate from \e,m 8 = 0) to 
\e,m s = d=l) and vice versa. P 9z i and P e ,i include spin pro¬ 
jections m s — ±1. (c) Energy levels and allowed transitions in 
NV center with radiatively-activated spin-mixing transitions. 
Km e and are the spin-mixing radiative excitation and 

decay rates. Spin-mixing of radiative origin [replacing K m in 
(b)] is described through K me and K mf rates that scale with 
the Purcell Factor (PF). 


be the random variables representing the number of pho¬ 
tons emitted during a measurement of | g,m s = 0) and 
| g,m s = ±1) respectively. Nq is the number of pho¬ 
tons emitted during a measurement of the NV spin state 
| g,m s = 0) and Ni is the number of photons emitted dur¬ 
ing a measurement of the NV spin state | g : m s = ±1). 
Since the number of photons emitted from the laser is 
Poisson distributed, we will assume that Nq and Ni are 
both Poisson distributed as well. Therefore the difference 
A = TVo — Ni is Skellam distributed m ■ Denoting no 
and ni to be the expected values of Nq and respec¬ 
tively, the expected value and variance of the difference, 
A, are 

E( A) = n 0 - ni, 

a 2 (A) = n 0 + ni. (1) 

The SNR of the Skellam distribution, defined as SNR = 
§[sy, is thus QU 

SNR = ^°~ ni . (2) 

Vno + ni 


The importance of the SNR stems also from its relation 
to NV state estimation. An (unexcited) NV state can be 
written as |Vv) = r\g,m s = 0) + (1 — r) | g : m s = ±1), 
where 0 < r < 1, is the fraction of ground state with 
spin projection m s = 0 (this definition could relate to a 
mixed-state or to a pure state, and, without loss of gener¬ 
ality, r is taken to be real for simplicity). In the Supple¬ 
mentary Material 03 we derive the maximum likelihood 
estimation of r and show that the minimization of the re¬ 
sulting estimation error is equivalent to maximizing the 
SNR in Eq. [2] 

In order to better understand the role of each of the 
system’s rates in the SNR, we present a theoretical anal¬ 
ysis of the NV measurement process. We rely on the 
following assumptions: 

1. The orbital mixing rate of the excited state is much 
larger than the decay rate [13] (this causes orbital 
averaging, essentially resulting in a single spin- 
triplet in the excited state). 

2. The radiative decay rate of the excited states is 
independent of the spin projection m- 

3. The excitation rate from both ground states is the 
same [T8] . 

4. The transition from the excited state with m s = 0 
to the singlet is negligible (since it is 4 orders of 
magnitude smaller than from m s = ±1) p~8l [T9] . 

5. Both ground state spin sublevels m s = —1 and 
m s = +1 have the same dynamics. 

6. The singlet state decays to | g,m s = 0) [IS] - 

7. The spin-mixing rate is the same between ±1 o 0. 

Given these assumptions, the rate equations that gov¬ 
ern the transitions in the NV center are: 

Pg,0 = ~K e Pg,0 4 " KfPe ,0 + N()Ps 

Pg, 1 = —KePg, 1 + KfP e ,l 
P s,0 = KePg, 0 — {Kf + 2iF m )P e? o + K m P e i 
Pe, i — K e Pg , i — (Kf + K s + K m )P e ^ i + 2iF m P e? o 
Ps = K s p eil - K 0 P S , (3) 

where P g ? o and P Q: i are the populations of the ground 
state with m s = 0 and m s = ±1 respectively, P e? o and 
P e? i are the populations of the excited state with mn s = 0 
and m s = ±1 respectively, and P s is the population of 
the singlet (Fig. [l]). K e is the excitation rate (spin- 
state preserving), Kf is the radiative decay rate (spin- 
state preserving), Kq is the decay rate from the singlet 
to \g,m s = 0), and K s is the decay rate from |e, m 3 = 1) 
(excited state with spin projection mn s = ±1) to the sin¬ 
glet. K m denotes small spin mixing transitions assumed 
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to be related to phononic coupling in the excited state 
pig. 0b)] m- The sum of the populations is normal¬ 
ized to be one, P g ? o + P Pj 1 + P e ,o + P e ,i + P s = 1. In the 
case of spin-mixing rates of radiative origin [Fig. [ljc)] 
the rate equations will be slightly different [Eq. (1) in 
m, with spin-mixing rates which scale with PF and K e . 
In the main part of the paper we assume non-radiative 
spin mixing, i.e. constant mixing rate, K m , between the 
excited states [Fig. 0b)]. 



FIG. 2. (a) SNR as a function of the pulse duration, T, for 
typical rates (see text). The optimal pulse duration is T — 
0.9911 [PL], at which the maximum value of SNR = 1.2516 
is reached, (b) SNR as a function of the excitation rate, K e , 
for typical rates and pulse duration T = 1[SL]. 


In the following analysis we used dimensionless vari¬ 
ables by normalizing the rates and durations by the life¬ 
time of the singlet state (~ 300 ns), which will be re¬ 
ferred to as the Singlet Lifetime (SL). Following this 
normalization the typical dimensionless NV rates are: 
Kf 0 = ^ « 23.077[irr] (the unmodified radiative de- 
cay rate), K s = (f^ -^) « 15.3846[^], K 0 = 1[^], 
K m = 0.0404iF/ o « 0.9323[^] [18]. We also take 
K e = Kf 0 as the typical excitation rate. 

For a given set of rates, the SNR, derived from the 
solution of the rate equations (Eq. [3|, peaks at an op¬ 
timal pulse duration, as is evident from Figure 0a) for 
the typical rates above. Moreover, for a given pulse du¬ 
ration the SNR increases as a function of the excita¬ 
tion rate, K e , and at large values of K e the SNR sat¬ 
urates asymptotically [Fig. |2jb)]. Both phenomena are 
in agreement with previously published experimental re¬ 
sults pH DE [20], thereby supporting our analysis and 
assumptions. For the typical NV rates the optimal pulse 
duration is T « 1 [SL] [Fig. [2ja)], which results in the 
maximum SNR = 1.2516. 

When adding a resonant optical structure to the sys¬ 
tem, the radiative decay rate, iLj, can also be con¬ 
trolled. The Purcell factor can be defined by the ratio 
of the modified decay rate, Kf, to the original rate, Kf 0 : 
PF = Kf/Kf 0 . Increasing the radiative decay rate in¬ 
creases both No and TVi, which have opposing effects on 
the SNR (since the contrast decreases). Figure 0a) de- 
picts the SNR as a function of the Purcell factor, PF, 
with typical NV rates and T = 1 [SL\. The plot clearly 




FIG. 3. The effect of Purcell enhancement on the SNR. (a) 
SNR as a function of PF for typical rates and T — 1. The 
SNR reaches a maximal value of SNR — 1.4867 at the op¬ 
timal Purcell factor PF « 3. The red dot shows the SNR 
for PF = 1. (b) SNR as a function of PF and T for typical 
rates. The SNR reaches a maximum value of SNR = 1.5951 
for optimal Purcell factor PF = 5.1903 and pulse duration 
T =s 1.703 [ST]. 


demonstrates that these conflicting effects give rise to 
a peak in the SNR that is achieved at an optimal PF. 
Specifically, the maximal SNR = 1.4867 is reached at 
the optimal Purcell factor PF « 3. Figure ^b) depicts 
the SNR as a function of both T and PF for typical NV 
rates, and demonstrates that optimizing the SNR over 
both PF and T results in a slightly higher SNR of 1.5951 
[compared to maximum SNR of 1.4867 with T = 1 [SL\ 
in Figure [3J a)], which is achieved for PF « 5.1903 and 
T = 1.703[ST]. This is a ^ 27% increase compare to the 
maximum SNR = 1.2516 with PF = 1 [Fig. |2ja)] . 




(a) (b) 

FIG. 4. Saturated SNR for non-radiative spin-mixing transi¬ 
tions. (a) The optimal pulse duration, T, as a function of the 
PF for K e approaching infinity, (b) Saturated SNR as a func¬ 
tion of the PF. The saturated value was calculated by taking 
the optimal T (for maximizing SNR) for each K e and calcu¬ 
lating the maximal SNR when K e goes to infinity. The red 
dot shows the saturated SNR with PF = 1. Increasing the 
PF and K e and adjusting T to its optimal value can increase 
the SNR significantly. 


We observe that changing the excitation rate, K ei or 
the Purcell factor alone does not increase the SNR consid¬ 
erably, as shown in figures [2] and [3j respectively. However, 
one might expect that optimizing T, while taking K e to 
infinity (SNR saturation value) as well as increasing the 
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PF, would take better advantage of the NV’s saturation 
behavior and thus increase the SNR significantly. For a 
fixed Purcell factor, both the optimal T (that maximizes 
the SNR) and the corresponding maximal SNR saturate 
when K e goes to infinity (see Fig. 1 in the Supplemen¬ 
tary Material HZD- This behavior allows us to investigate 
the effect of the Purcell factor on the asymptotic value of 
the maximal SNR (which will be referred to as the satu¬ 
rated SNR), and the corresponding asymptotic value of 
T (which will be referred to as the optimal T). We note 
that choosing K e = 2 Kf = 2 PF * Kf 0 approaches the 
saturation SNR values to within 15%. 

Fig. [4] depicts the saturated SNR as a function of PF 
for optimal T, assuming non-radiative spin-mixing rates 
that are not affected by the change in PF. The results 
clearly show a significant increase in the saturated SNR 
when increasing the PF, with no observable saturation. 
This demonstrates that under the assumptions that the 
spin mixing terms are of phononic origin (i.e. do not 
scale with the optical rates), the SNR can be improved 
significantly by increasing PF. This is the main result of 
this work. In particular, by reaching a moderate PF of 
only 4 the SNR can be doubled (compare to PF = 1). 
Such a low PF should be easily achieved using broadband 
plasmonic nano-antennas for example, emphasizing the 
feasibility of such an SNR enhancement scheme. 




(a) (b) 

FIG. 5. Saturated SNR for radiative spin-mixing transitions, 
the mixing rates were set to be Km f = 0.02 Kf « 0.4615PF 
and Km e — 0.02 K e . (a) The optimal pulse duration, T, as a 
function of the PF for K e going to infinity, (b) Saturated SNR 
as a function of PF. The saturated value was calculated by 
maximizing the SNR over T as K e approaches infinity. The 
red dot marks the saturated SNR with PF = 1. Optimizing 
Kf , K e and T results in a maximal SNR value at the optimal 
PF « 2. 


We now analyze the opposite scenario, in which the 
spin-mixing terms are of radiative origin [Fig. [ljc)] . 
Figure [5] shows that in this case the saturated SNR is 
bounded, and reaches its maximum at PF ~ 2. Optimiz¬ 
ing the saturated SNR results only in a ^ 6% improve¬ 
ment over the maximum SNR with the original decay 
rate (PF = 1). The fact that the spin-mixing process 
is optically activated, and thus scales with K e and PF, 
strongly suppresses the SNR enhancement (since the spin 


contrast is limited). 

We note that the measurements proposed in this work 
are highly sensitive to the relatively small spin-mixing 
terms. We find a significant (unsaturated) increase in the 
SNR assuming constant spin-mixing rates (Fig. [4]), while 
for radiative mixing rates (which scale with the Purcell 
factor) the results indicate only a moderate increase in 
the spin-readout SNR (Fig. [ 5 ]). This strong dependence 
on the unknown origin of the spin mixing could provide 
a useful approach for studying the underlying physical 
processes related to this mixing at room temperature [20]. 

Changing the different rates of the system may also 
affect the optical initialization of the NV to | g,m s = 0). 
However, increasing the PF has only a minor effect on the 
ability to initialize the NV center. Increasing the PF to 
10 causes a polarization decrease of only ^ 1% after ini¬ 
tialization in the case of non-radiative spin-mixing, and 
a ~ 6% decrease for radiative spin-mixing rates. 

In conclusion, we have analyzed the spin-readout SNR 
for NV centers in diamond, in the presence of optical 
couplers affecting their fluorescence rate through Purcell 
enhancement. We have constructed a relevant measure 
of the readout SNR, taking into account two competing 
processes - elevated signal yet reduced signal contrast - 
with increasing radiative decay rates. We find that the 
SNR can be improved by concurrently increasing both 
the decay and excitation rates, such that the NV satura¬ 
tion behavior is fully exploited. We find a significant and 
unsaturated improvement in the SNR as a function of the 
Purcell factor, assuming that spin-mixing rates are con¬ 
stant (unaffected by the Purcell enhancement). We note 
that this large improvement in readout SNR does not 
take into account potential improvements in collection 
efficiency that could result from the use of such optical 
couplers (e.g. mm), and could have a major impact 
on a variety of applications, potentially leading to single¬ 
shot spin readout. We also analyzed the system assuming 
that spin-mixing is radiative and scales with the Purcell 
factor. In this case the SNR is severely limited and sat¬ 
urates at ~ 6% above its natural value (with PF = 1). 
The strong dependence of our results on the origin of the 
spin-mixing terms suggests this approach as useful for 
studying the physics dominating such processes. 

While preparing this manuscript we have become 
aware that related work is described in the thesis [22 . 
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